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We report on the control of the faceting of crystal surfaces by means of surface electromigration. 
When electromigration reinforces the faceting instabihty, we find perpetual coarsening with a wave- 
length increasing as t^'"^ . For strongly stabilizing electromigration, the surface is stable. For weakly 
stabilizing electromigration, a cellular pattern is obtained, with a nonlinearly selected wavelength. 
The selection mechanism is not caused by an instability of steady-states, as suggested by previ- 
ous works in the literature. Instead, the dynamics is found to exhibit coarsening before reaching a 
continuous family of stable non-equilibrium steady-states. 
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In non-equilibrium conditions, crystal surfaces undergo 
various instabilities leading to micro or nanostructures. 
The ultimate faith of these topographical structures is 
governed by the nonlinear dynanrics of the surface. Un- 
derstanding the nonlinear processes at play is therefore 
important for many applications where one wishes to de- 
sign structures with a given wavelength. However, the 
complexity of nonlinear wavelength selection mechanisms 
have hindered their control, and strategies to achieve pre- 
dictive design arc still poorly understood. A possible 
strategy is to control existing instabilities by means of an 
external field. Here we explore the control of the faceting 
instability with electromigration. The faceting instabil- 
ity is the decomposition of a surface with a given aver- 
age orientation into neighboring facetted orientations [l|- 
|5|. This instability is of thermodynamic origin, and is 
driven by the reduction of the surface free energy. The 
presence of an electric current leads to a surface elec- 
tromigration mass flux Je- This mass flux can stabi- 
lize or destabilize crystalline surfaces depending on its 
orientation-dependence [6|-|8|]. Many studies have been 
devoted to these two instabilities. The combination of 
electromigration and anisotropy is already known to lead 
to nontrivial dynamics of voids, with the appearance of 
nontrivial orientations and oscillatory dynamics [9|. In 
this letter, we show that a faceting instability can be 
controlled by electromigration. In the presence of desta- 
bilizing electromigration, the faceting instability is found 
to be reinforced, and perpetual coarsening with a wave- 
length increasing as t^/^ is found. Strongly stabilizing 
electromigration supersedes the faceting instability and 
stabilizes the surface. This result is similar to the sta- 
bilization of the clastic stress-induced Grinfeld instabil- 
ity by electromigration discussed in Ref. [10| . Our most 
striking result appears in the presence of weakly stabi- 
lizing electromigration: we provide converging analytical 
and numerical evidences of the existence of a continuous 
branch of stable periodic steady-states. Since the dynam- 
ics exhibits a Lyapunov functional [i.e. an effective non- 
equilibrium energy functional which is decreased mono- 



tonically during the dynamics), these steady-states can 
be characterized as non- equilibrium meta-stable states. 
Furthermore, when starting from random initial condi- 
tions, a wavelength larger than the wavelength emerging 
from the linear instability is selected. However, the non- 
linear wavelength selection mechanism is different from 
the known scenarios of interrupted coarsening discussed 
in Refs. llMl3l|. or of secondary instabilities such as the 
Eckhaus instability |14{ . which all emerge from instabili- 
ties of steady-states. 

Model - Combining previous works on faceting [J, |5| 
and on clectronrigration- induced instabilities |6|-|8|, we 
write down a phcnomenological nonlinear equation which 
governs surface dynamics. We use a one-dimensional 
model for the crystal surface height profile h{x,t) and 
slope (j){x,t) = dxh{x,t). The average orientation = 
of the surface is assunred to be unstable, and to decom- 
pose into facets of slopes (f> = ±1. For the sake of sim- 
plicity, we consider a Ginzburg-Landau-like orientation 
dependent energy: 



J^[0] ^ ■^ dx 



-1 + ^ + ^(4*)^ 



(1) 



where 7 is a typical surface energy scale. The last term 
on the r.h.s. accounts for a curvature energy cost, which 
regularizes facets [J, l5| at the small lengthscale e, and 
allows one to write down local dynamics. From Eq.(IT|), 
the local chemical potential is derived as 



/i = il 



5h 



l^dx [e'5. 



(2) 



where Vt is the atomic area. We consider surface-diffusion 
limited dynamics, so that chemical potential gradients 
induce a surface mass flux J^ = —DLdxli/kBT, where 
ksT is the thermal energy. In addition, we assume that 
an external electric current is applied to the crystal, lead- 
ing to an orientation-dependent electromigration surface 
mass flux Je{4'). Expanding the electromigration cur- 
rent for small slopes, we write Jsifj)) ~ Je{0) + 4>J'e{0)- 
The slope-dependence of Je uiay either be caused by the 



anisotropy of surface diffusion [6|, |l5|, or of the migra- 
tion force [Ig. As shown in Fig[T{a), electromigration 
is destabihzing when J^(0) > 0, and stabihzing when 
J'e{0) < 0. This instabihty is the subject of a large hter- 
ature, and has been the basis of the interpretation of the 
instabihties observed on semiconductors 17[ and on met- 
als [18[. Mass conservation reads dth = —fldx[Jfj, + Je], 
and leads to an evolution equation for the surface height: 



dth 



nj 



kbT 
J'EWdxxh 



C'xx [^ Oxxxx'^ \ Oxx*^ (^xxpx'^) 



(3) 



where we have assumed for simplicity that D^ docs not 
depend on the orientation. Normalizing x with e, t with 
kBTe^ /VL^^iDl, and defining the dimcnsionlcss param- 
eter j'e = JE{0)kBTe^ /QjDl, we may write a one- 
parameter equation for the slope: 



dt<l> = d^x^^ [dxx<l> + f/* - 0^] - Je9xx(I) 



(4) 
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The relaxation part of Eq. ^ , first derived in Ref . 
expected to give rise to p erpetual logarithmic coarsening 
as discussed in Ref. [l9[. The electromigration part has 
been derived and used in many papers, see Ref. [Sl for a 
review. When electromigration is destabilizing (j^ > 0), 
powerlaw coarsening is expected. Here, we claim that the 
combination of the two physical processes, faceting and 
electromigration, gives rise to novel dynamics. 

Linear stability analysis ~ Let us start the study of 
Eq.((3]) by means of a linear stability analysis. The growth 
rate iuj of a Fourier mode ~ Qii^t+qx ^^ wavelength A = 
27r/g reads: 



ioj = -g^ + q^+ j'eQ^ 



(5) 



For j'e = 0, all wavelengths larger than 27r are unstable. 
For Je > the range of unstable wavelengths is still in- 
finite, but extends to smaller and smaller wavelengths: 
A > A+ = 27r/[l/2+ (1/4 + j^)i/2]i/2. when electro- 
migration is weakly stabilizing > j^ > —1/4, there 
is a finite range of unstable wavelengths A_ < A < A+ 
where A_ = 27r/[l/2 - (1/4 + j^)!/^]^/^^ Finally, for 
strongly stabilizing electromigration, when j'^ < —1/4, 
the surface is linearly stable. When the surface is unsta- 
ble, the wavelength A„ = 2n3^/^/[l + (1 + 3j^)i/2]i/2 
of the fastest growing mode is expected to emerge from 
random initial conditions (where all modes are present). 
Lyapunov functional - The linear stability analysis 
predicts the existence or not of an instability from small 
initial perturbations. However, the amplitude of the un- 
stable modes increases exponentially fast, and the un- 
stable front enters into the nonlinear regime where the 
(j)^ term in Eq.Q is not negligible anymore. While the 
analytical study of the nonlinear regime is in general del- 
icate, the existence of an effective energy functional C, 
sometimes called a Lyapunov functional, usually greatly 
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FIG. 1: Mechanism of the instability and dispersion relation: 
(a) Schematics of the faceting and electromigration-induced 
instabilities, (b) Growth rate ioj of Fourier modes from linear 
stability analysis, see Eq. (O. 



simplifies the analysis. Such a functional is defined as 
a quantity which decreases monotonously during the dy- 
namics. The evolution equation Q actually exhibits a 
Lyapunov functional (in normalized coordinates): 



C[li] = T[dxh] 



Je 



dx ti{x) 



(6) 



Indeed, one may write dth = dxx[SC/Sh], leading to: 



dtC 



dx[dx{5C/5h)f <0. 



(7) 



Steady-state branches - The space of all possible con- 
figurations for the surface profile is very large, and it is 
impossible to calculate the value of C for all configura- 
tions. As a consequence, we would like to reduce this 
space to a relevant subset of shapes which is simpler 
to explore. A powerful approach along this line is to 
study the steady-states </>o(a^), which are the solutions of 
dt(j)o = 0. Using Eq.©, one has dx{SC/Sh) = 0, leading 
to 



dxx [dxx(t)Q + (1)0- (1)1] - j'E(f>a = 0. 



(8) 



The central idea motivating the study of steady-states 
is the existence of a separation of timescales, where the 
shape first relaxes rapidly towards periodic steady-states. 
Then, these periodic solutions may exhibit an instability 
occurring at longer timescales, leading for example to a 
coarsening process [l3| , or to the Eckhaus instability [iJ] . 
The steady-states profiles are obtained numerically, 
starting with a small amplitude sinusoidal perturbation 
~ sin(27ra:;/A) with one period in a box of width A|23l|. 



Following the linear stability analysis, the perturbation 
may grow or decay. When it grows, it reaches a finite 
amplitude steady-state where it stops. As shown in FigH] 
for j'^ > 0, the steady-state exhibits a monotonously in- 
creasing amplitude ^ as a function of the wavelength A. 
Asymptotically for large A and A, we expect the domi- 
nant term among the linear terms to be the one with the 
fewest derivatives, i.e. j'e4'0i ^-iid this term has to bal- 
ance the nonlinear term d^xfj^o^ leading to ^ ~ (Jb)^^^'^- 
The numerical determination of the steady-state branch 
actually indicates that: A w 0.16(j^)^/^A. 

When —1/4 < j'^ < 0, we obtain a bell shape branch 
connecting the two marginal points (A — X-,A — 0), 
and (A = A-)-, A = 0), as shown in Fig|3|a). In addition, 
using suitable initial conditions, we find another steady- 
state branch, which emerges from the primary bell-shape 
branch. The typical steady-state profiles are shown on 
FigEId). Other branches could exist, and we have not 
tried to provide a complete analysis of all possible steady- 
state branches. However, and as discussed in the follow- 
ing, the study of the main bell-shape branch seems to 
be sufficient to account for the main features of the full 
dynamics starting from small random perturbations. 

The numerical calculation of C along the steady-state 
branches shows that C decreases monotonously with A 
for j'^ > 0, while it exhibits a minimum for A = A^ when 
— 1/4 < j'^ < 0. At this point, it is tempting to speculate 
that the dynamics will simply follow the gradient of C 
along the steady-state branches [2J], leading to infinite 
coarsening when j'^ > and interrupted coarsening at 
A = A£ for —1/4 < j'^ < 0. We shall see in the following 
that some of these speculations are actually wrong. 

Steady-state stability - In order to analyze the dy- 
namics more carefully, we shall investigate the stabil- 
ity of steady-states with respect to small perturbations 
4>i{x) = (j>{x) — (l)o{x) ^ (poix). The slope variation (/>!, 
leads to the following variation of C: 



C^ = ^Jdx [(/.2(302 _ 1) + (a^0i)2 _ f^hl] 



(9) 



where dxhi = 0i. Physically relevant steady-states must 
be stable with respect to perturbations with wavelengths 
smaller than their periodicity, and in general, we expect 
that the most dangerous modes are long- wavelength per- 
turbations, as pointed out e.g. in Ref. [l3[. In the long 
wavelength limit, where the perturbation wavelength is 
much larger than the periodicity of the steady-state 4>o, 
one may simply replace (f>Q (x) by its average over one pe- 
riod (0o)- Thus, since it has constant coefficients, Eq.® 
is now diagonal in Fourier space, and this suggests a sim- 
ple stability criterion: 



(3(0^) - l)ql 



it 



Je > 0. 



(10) 




FIG. 2: Perpetual power-law coarsening for j'^ > 0; (a) The 
black line indicates the steady state amplitude A versus wave- 
length A of periodic steady-states for j'^ = 0.45. Crosses 
(x): full simulations from small random initial conditions. 
(b) Rescaled Amplitude A/j'^ and wavelength A/(j^)'^" as 
function of time for starting from small random initial con- 
ditions. The solid lines are the powerlaws discussed in the 
text. 



tion j'^ > 0. This is in agreement with the above specu- 
lation of perpetual coarsening. In addition, we may also 
gain information about the coarsening exponent. Indeed, 
at long wavelengths we have £i « — / dx j^/if /2, and as 
a consequence one has dthi ~ —j'^dxxhi. This relation 
provides a link between lengthscales x, assumed to be 
~ A, and timescales t, leading to A '-^ [j'eY^'^^^^'^- Using 
the previously derived linear relation between A and A^ 
we also find A ^ j^t^/^. 

When -1/4 < j^ < 0, from Eq.dTU]), the steady-states 
should be stable for all qi when 



^i)>^[i-2{-jEr^'] 



(11) 



As a consequence, long wavelength modes qi — >■ are 
always unstable in the case of destabilizing electromigra- 



This criterion indicates that the upper part of the branch 
between the green stars in Fig|3lja,b) should be stable, 
while the lower parts close to A± should be unstable [25 1 . 
Such a result suggests that the system is stuck once the 
dynamics hits the stable part of the steady-state branch, 
and as a consequence, it cannot evolve towards the min- 
imum of £ at A = A£ . 

Full dynamics - These results are confirmed by the 
full numerical solution of Eq. Q starting from small ran- 
dom perturbations of a flat state in a system of size 
L = 500. The results are shown in FigHl First, and 
as expected, we find stable dynamics for j'^ < —1/4. 

For destabilizing electromigration j'^ > 0, perpet- 
ual coarsening is found, as shown on FigUJ^a), and af- 
ter a transient related to the linear instability, the dy- 
namics follows the steady-state branch in the (A, A) 
plane, as shown in Figl2ja). We also confirm the scal- 
ing laws given above in Figl2jb), and find the prefac- 
tors: A « lA{j'j^)^/H^/^, and A « 0.22 j'j^t^/^. The ra- 
tio of these prefactors is in perfect agreement with the 
asymptotic behavior of the steady-states discussed above 
A«0.16(j^)i/2A. 
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FIG. 3: Nonlinear wavelength selection for —1/4 < j'^ < 0: 
(a) The (black) solid and (brown) dashed lines indicate the 
steady state amplitude A versus wavelength A for j'^ = —0.15. 
Red line: full dynamics from small random initial conditions 
in an extended system. Orange lines: dynamics starting from 
slightly perturbed periodic steady-state in an extended sys- 
tem. Green stars indicate theoretical limit of stability of the 
steady-state branch, (b) The shaded region corresponds to 
the linearly unstable region. Green stars: limits of steady- 
state stability, (c) Lyapunov functional density C/X evalu- 
ated from the steady state profiles, (d) Steady state profiles 
numerically obtained for the two branches at j'^ = —0.15. 



For weakly stabilizing electromigration —1/4 < j'^ < 
0, the dynamics starts with a rapid increase of the am- 
plitude due to the linear instability. Then, the average 
wavelength increases before the dynamics hits the steady- 
state branch at A = A*, as shown in Fi^JlJa). When the 
dynamics reaches the steady-state branch, it stops, as 
expected from the above prediction of steady-state sta- 
bility. As a consequence, the system never reaches the 
minimum of C along the steady-state branch. 

In order to check further the stability of the steady- 
state branch, we have performed simulations in boxes of 
width around 80 periods starting with periodic steady- 
states with small perturbations. These simulations con- 
firm the stability of the upper part of the branch, and 
the instability of the lower part, in quantitative agree- 
ment with the condition ([TT|) . As seen in Figl^lja), the 
instability does not lead to a trajectory of the system 
along the steady-state branch in the {X,A) plane. In- 
stead, the trajectory escapes from the branch and returns 
to it, stopping in the same region as the dynamics from 
flat initial conditions. Fig|3l^b) summarizes the evolution 
of the different lengthscales as a function of j^. 

Conclusions and perspectives - In summary, we have 
studied the control of the faceting instability by means of 
electromigration. For strong stabilizing electromigration. 
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FIG. 4: Full dynamics and coarsening: Left: surface height 
h{x, t) as a function of space for different times in the three 
different regimes: (a) j'^ = 0.45, (b) j'^ = —0.15, (c) j^ = 
—0.3. Right: Corresponding spatiotemporal portrait of the 
extrema of the height profile. 



the surface is stable. Under weakly stabilizing electromi- 
gration, the surface exhibits a periodic cellular structure 
with a nonlinearly selected wavelength. When electro- 
migration is destabilizing, perpetual coarsening is found 
with a coarsening exponent 1/2. 

In the case of weakly stabilizing electromigration, we 
find a continuous family of stable periodic steady-states. 
These states can be interpreted as non-equilibrium 
metastable states, where metastability is here defined 
as local stability with respect to a non-equilibrium Lya- 
punov functional playing the role of an effective energy. 

In the literature, several cases of nonlinear wavelength 
selection on crystal surfaces have been interpreted as in 
terrupted coarsening, such as for mound growth 
atomic step meandering [1^, and ion sputtering 



111, 
2C|. 



However, the nonlinear wavelength selection scenario pre- 
sented in this Letter does not correspond to interrupted 
coarsening, as defined e.g. in Ref. [13[, or to another in- 
stability of steady-states, such as in the Eckhaus instabil- 
ity UM- Indeed, we only observe a small amount of coars- 
ening before the dynamics hits the steady-state branch. 
Once the system has reached a steady-state, it is stable 
and the evolution stops. 

These results may provide hints to understand other 
nonlinear wavelength selection scenarios obtained in the 
literature, such as in the combination of growth and 
faceting as discussed in Ref. [2l|, [2^ (where nonlinear 
wavelength selection occurs between a coarsening regime, 
and a chaotic regime) . More generally, we hope that our 
work will provide milestones towards the novel methods 
to control the size of nano-structures emerging from sur- 
face morphological instabilities. 

KM was supported by the Marie Curie FP7-PE0PLE- 
2009-IEF program. 



[is; 



[1] 

[2] 

[3] 

[4] 
[5] 
[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 



A. Chame, S. Rousset, H. Bonzel, and J. Villain, Bulg. [14 

Chem. Commun. 29, 398 (1996). 

Y. Kim, M. H. Jo, T. C. Kim, C. W. Yang, J. W. [15 

Kim, J. S. Hwang, D. Y. Noh, N. D. Kim, and [15^ 

J. W. Chung, Phys. Rev. Lett. 102, 156103 (2009), URL 

http : //link . aps . org/doi/10 . 1103/PhysRevLett . 102 . 156103| 



[17; 
[is; 

[19 



H.-C. Jeong and J. D. Weeks, Phys. 
Rev. Lett. 75, 4456 (1995), URL 

[http : //link . aps . org/doi/10 . 1 103/PliysRevLett . 75 ■ 4456 [ 
A. Di Carlo, M. Gurtin, and P. Podio-Guidugh, SIAM J. 
Appl. Math. 52, 1111 (1992). 

F. Liu and H. Metiu, Phys. Rev. B 48, 5808 (1993), URL 
,littp : //link . aps . org/doi/10 . 1103/PhysRevB . 48 . 5808 , 
J. Krug and H. T. Dobbs, Phys. 
Rev. Lett. 73, 1947 (1994), URL [20 
|http : //link . aps . org/doi/10 . 1103/PhysRevLett . 73 . 1947[ 

J. Chang, O. Pierre-Louis, and C. Misbah, 
Phys. Rev. Lett. 96, 195901 (2006), URL 
[http : //link. aps . org/doi/10 . 1103/PhysRevLett . 96 . 19590 Ij i 
C. Misbah, O. Pierre-Louis, and Y. Saito, 
Rev. Mod. Phys. 82, 981 (2010), URL 
'http : //link . aps . org/doi/10 . 1103/RevModPhys . 82 . 981 ' [22 
P. Kuhn, J. Krug, F. Hausscr, and A. Voigt, 
Phys. Rev. Lett. 94, 166105 (2005), URL 
,littp : //link . aps . org/doi/10 . 1103/PliysRevLett . 94 . 166105 1 
V. Tomar, M. R. Gungor, and D. Maroudas, r23 
Phys. Rev. B 81, 054111 (2010), URL 34^ 

|http : //link . aps . org/doi/10 . 1103/PhysRevB . 81 . 0541111 
P. Politi and J. Villain, Phys. 

Rev. B 54, 5114 (1996), URL 

|http : //link . aps . org/doi/10 . 1103/PhysRevB . 54 . 5114] 

G. Danker, O. Pierre-Louis, K. Kassner, and 
C. Misba h, Phys. Rev. E 68, 020601 (2003), URL 

|http : //link. aps . org/doi/10 ■ 1103/PhysRevE . 68 ■ 020601 [ 



[25] 



P. Politi and C. Misbah, Phys. 

Rev. Lett. 92, 090601 (2004), URL 

|http : //link . aps . org/doi/10 ■ 1103/PhysRevLett . 92 . 090601 | 

W. Eckaus, Studies in Nonlinear Stability Theory 

(Springer, Berlin, 1965). 

S. Stoyanov, Jpn. J. Appl. Phys. 29, L659 (1990). 

O. Pierre-Louis, Phys. Rev. 

Lett. 96, 135901 (2006), URL 

|http : //link . aps . org/doi/10 . 1103/PhysRevLett . 96 . 13590T] 
A. Latyshev, A. Aseev, A. Krasilnikov, and S. Stenin, 

Surface Science 213, 157 (1989), ISSN 0039-6028, URL 

|http: //www. sciencedirect . com/science/article/pii/0039602889 

R. P. Johnson, Phys. Rev. 54, 459 (1938), URL 

http : //link . aps . org/doi/10 . 1 103/PhysRev . 54 . 459 

F. Haufier and A. Voigt, Phys. 

Rev. E 79, 011115 (2009), URL 

|http : //link . aps . org/doi/10 ■ 1 103/PIiysRevE . 79 . OliTlS] 

J. Mufioz Garcia, R. Gago, L. Vazquez, 

J. A. Sanchez-Garci'a, and R. Cuerno, 

Phys. Rev. L ett. 104, 026101 (2010), URL 

|http: // link. aps . org/doi/10 . 1103/PhysRevLett . 104 .026101] 

A. A. Golovin, A. A. Nepomnyashchy, S. H. Davis, and 

M. A. Zaks, Phys. Rev. Lett. 86, 1550 (2001), URL 

|http:/ /link . aps . org/doi/10 . 1103/PhysRevLett . 86 . 1550[ 

T. V. Savina, A. A. Golovin, S. H. Davis, 

A. A. Nepomnyashchy, and P. W. Voorhees, 

Phys. Rev. E 67, 021606 (2003), URL 

|http : //link . aps . org/doi/10 . 1 10 3/PhysRevE . 67 ■ 021606 | 

We use an finite diflcrcncc explicit scheme. 

More precisely, since A varies, we calculate the density 

C/X. 

Due to numerical difficulties in the calculation of the 

steady-state branch for A close to A+ when > j^; > 

—0.1, the upper stability limit in Fig. 3(b) was obtained 

from an interpolation for j'^ = —0.1, —0.05. 



